Effect of Coulomb interactions on the physical observables of graphene 
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We give an update of the situation concerning the effect of electron-electron interactions on the 
physics of the neutral graphene system at low energies. We revise old renormalization group results 
and the use of the 1 /N expansion to address the questions of the possible opening of a low energy 
gap, and the magnitude of the graphene fine structure constant. We emphasize the role of the 
Fermi velocity as the only free parameter determining the transport and electronic properties of 
the graphene system and revise its renormalization by Coulomb interactions on the light of recent 
experimental evidences. 

PACS numbers: 

I. INTRODUCTION 

Graphene is considered to be a bridge between quantum field theory (QFT) and condensed matter physics due to 
the low energy description of the elementary excitations as massless Dirac fermions in two spacial dimensions. This 
description realizes some important toy models used to study quark confinement and chiral symmetry breaking in 
QFT, particularly Quantum Electrodynamics in two spacial dimensions QED(3) [H,i. On the condensed matter 
point the QFT modeling of graphene raises two important questions: under a fundamental point of view it is not 
guaranteed that the neutral graphene system with Coulomb interaction is a Fermi liquid [3]. On a more practical 
side, the QFT model has logarithmic divergences that need to be renormalized what introduces some subtleties in 
fixing the value of the graphene fine structure constant affecting the experiments. 

Electron-electron interactions were considered to be very small or even absent in the first experiments 0413 but 
the recent advances in the synthesis and the big improvement in the quality of the samples have changed this view. 
The observation of the fractional Quantum Hall effect [TH, [l4[ and the experiments probing the zero Landau level 
in high magnetic fields fl5{ have renewed the interest on the role of Coulomb interactions on the intrinsic properties 
of the system. Since the role of many body corrections to the physics of graphene has been extensively studied in 
the literature (for a recent review see [lH and the references inside), we will here pick up these aspects that in our 
opinion are more interesting or less well understood. Here we just note that, since all other interactions are irrelevant 
at sufficiently low energies [I?], HH > the Coulomb interaction is the only feature affecting the physical properties of 
the clean suspended samples. 

The low energy effective model for graphene being a QFT model suffers from one typical problem: perturbation 
theory calculations of any observable quantity are plagued by infinities or spurious cutoff dependences. The standard 
QFT renormalization program allows us to address two very important questions in the physics of graphene: what is 
the strength of the Coulomb interaction at low energies - i. e. what is the actual value of the graphene fine structure 
constant -, and will the interactions lead to a band gap opening? We summarize the contents of this review by giving 
answers to these questions, the second first: A gap will not open in graphene at low energies for moderate values of 
the graphene fine structure constant ac- And: the infrared value of qq is smaller than the assumed value based on 
a Fermi velocity of the order of c/300. 

Note however that interactions are very much enhanced in the presence of strong magnetic fields. In this situation 
a gap can open at the neutrality point [19J. This article describes the physical low energy properties of neutral, clean 
graphene at zero magnetic field. 



II. GRAPHENE VERSUS QED 
A. The continuum model 



One of the facts that has allowed the fast development in the understanding of the electronic properties of graphene 
is that there is a well established model Hamiltonian universally accepted by the community. This is another similarity 
with the physics of elementary particles. The standard non-interacting model for the electronic excitations around a 
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single Fermi point in graphene (in units h = 1) is given by 

H = vf /d 2 r^(r)ya^(r) , (1) 



where i = 1,2, 4>{r) — V ,+ (r)7°, and the gamma matrices can be chosen as 7^ = 02,7,, = — o"i,7 = 0-3 • c,; are 
the Pauli matrices and vp is the Fermi velocity, the only free parameter in the Hamiltonian. When promoting the 
Hamiltonian to a QFT action we have the electron wave function as an additional parameter defining the theory. 
And the electron-electron interactions are described by the Hamiltonian 

nint ^ e , f dW tim^ipm. (2) 



It is immediately seen that, unlike what happens in the usual two dimensional electron gas, the ratio of the Coulomb 
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interactions and the kinetic energy in this system is a constant independent of the density and given by g ~ ^ often 
taken as the graphene fine structure constant. 

Since the model Hamiltonian is the relativistic Dirac Hamiltonian it saves a lot of efforts to adopt a covariant 
formulation and profit of the results already known in QFT. The fact that the numerical value of Fermi velocity does 
not equal the speed of light is not important, the model is still relativistic with a different light cone. What is very 
different and interesting under a QFT point of view is the fact that the velocity defining the light cone is renormalized 
by the interactions. 

Following the quantum field theory nature of the model, we replace the instantaneous Coulomb interaction of eq. 
@ by a local gauge interaction through a minimal coupling: 

Ant =9 J #xdtj' 1 {x,t)An(z,t) , 

where the electron current is defined as 

f = (*7 *,« J F*7 i *) . 
With this substitution we have formally an effective Lagrangean 

C = [ d 2 vdt^ 7^ + ieA^ , (3) 



very similar to that of massless QED in two spatial dimensions. What makes the graphene model different from QED 
in any number of dimensions is the gauge boson propagator. As discussed at length in [13, 0, H3, a gauge field 
propagator in QFT has a 1/k 2 dependence while the effective Coulomb interaction in graphene propagates as 1/k. 
This is due to the fact that only the charges are confined to live in the two dimensional surface defined by the sample. 
The photons propagate in the three dimensional space. This is a crucial aspect that makes the model interesting and 
behaving more like the scale invariant massless QED(3+1) rather than the super-renormalizable QED(2+1). The 
renormalization of the Fermi velocity is a consequence of this fact. 



B. Perturbative renormalization 



In standard QED (4) there are three one loop diagrams that have logarithmic singularities (see Fig. [1} and three 
free parameters in the theory (the coupling constant and the electron and photon wave functions). The theory is 
strictly renormalizable in the sense that all divergences at any order in perturbation theory can be cured by a proper 
redefinition of the parameters. Subsequently the physical values of these parameters are fixed by experimental data 
(renormalization conditions). On the other hand, QED(3) is a super-renormalizable theory meaning that the coupling 
constant has dimensions of energy what improves the convergence of the perturbative series. It has less divergences 
and these can be properly renormalized. In fact, massless QED (2+1) is ultraviolet finite although it has infrared 



divergences [2l|. Graphene sits in between QED(3) and QED(4). Of the three graphs shown in Fig. [T] only the 



electron self-energy diverges and in the case of considering a static photon propagator only the spatial part has a 
logarithmic divergence. In the static model, the electron and photon propagators in momentum space are given by 



(4) 



b) 





FIG. 1: Primitively divergent Feynman graphs in QED(4). Fermion self-energy (a), vertex (b) and photon self-energy (c) 



n (k) = ^. (5) 

The renormalization functions associated to the electron self-energy are defined as: 

Go 1 - E(*°,k) = Z; 1/2 (fc°,k)[fc° 7 - Z v (k°,k)v 7 ■ k], (6) 

The extra parameter v appearing in ((4]) in the graphene case has an associated renormalization function Z v which is 
a new feature characteristic of graphene. In the Lorentz invariant massless QED in any dimensions, only the electron 
wave function renormalization is associated to the the electron self energy. The computation of the graph in Fig. [T] 
(a) gives 

£«(k)~|« 7 .klog^, (7) 

where g is the coupling constant to be discussed later, and A is a high energy cutoff. The logarithmic divergence can 
be fixed by redefining the fermion velocity v which becomes, after finishing the renormalization procedure, energy- 
dependent. The result is by now well known: 

^ = l--f ^ g (E/E ), (8) 

VF(Eq) I07T 

which relates the value of the Fermi velocity at an energy E, with that at a reference energy Eq, assuming that 
the two energies are sufficiently closed. Notice that Eg is not a cutoff of the order of the bandwidth but any 
energy where the Fermi velocity is determined by experiments. The corresponding diagram in massless QED(3) 
does not have a logarithmic divergence because the photon propagator has an extra inverse power of momentum. In 
QED (4) this diagram induces a wave function renormalization. It is interesting to note that the electron charge in 
the graphene model is not renormalized (after fixing the velocity divergences, the photon self-energy is finite at all 
orders in perturbation theory) but the velocity renormalization induces a renormalization of the graphene structure 
constant. Since according to eq. flSJ) the Fermi velocity increases at decreasing energies, the effective coupling constant 
g = e 2 /4ttvf decreases in the infrared making perturbation theory more accurate. 

III. THE COUPLING CONSTANT OF COULOMB INTERACTION AND PERTURBATION THEORY 

As mentioned in the introduction perturbation theory calculations of any observable quantity are plagued by 
infinities. Perturbative renormalization |22j is a well prescribed mechanism to get rid of the infinities and define 
a sensible predictive theory. In the graphene system we encounter an additional problem: The bare value of the 
graphene fine structure constant 

«G = ~7~~ ' (9) 

is not small. Plugging-in the value of the bare electron charge e and the Fermi velocity measured in nanotube 
experiments [23[ vf ~ c/300, etc is estimated to be of the order of ag ~ 2.3 — 2.5. 

A. Screening in graphene 



The estimate of a described above precludes the use of perturbation theory. This problem is not so severe for 
samples on a substrate, where the substrate dielectric constant reduces the effective coupling constant: 
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where eg includes intrinsic contributions and effects due to the environment in which graphene is immersed. A simple 
estimate gives Eq = (£ a ir + £subs)/2 ~ 2 — 3 for typical substrates. An interesting topic is the intrinsic screening due 
to the graphene layer itself. Internal excitations in three dimensional insulators lead to a finite dielectric constant, 
which reduces the effective charge of the mobile carriers. By extrapolating measurements of the excitation spectra 
in graphite, it has been proposed that similar effects lead to a large intrinsic dielectric constant in single layer 
graphene [24}. eg ~ 13. 

It is interesting to note that internal excitations in insulators, including electron-hole pairs, can be described as 
polarizable dipoles, which induce a field £md(q) oc |q| 2 £ e2 . t (q), where q is the momentum of the external field. In a 
two dimensional system dipoles induce a field which tends to zero at small momenta, so that they do not contribute 
to the dielectric constant. The only internal screening processes in graphene which can modify the dielectric constant 
are the excitations near the Dirac energy. Their one loop contribution to the polarizability is finite and it can be 
calculated analytically 

e 2 k 2 

n(k,a;)=i- (11) 

8 yJv 2 p\L 2 — LO 2 

This one loop result is independent on the nature of the interaction since only electron propagators appear in the 
calculation. The resulting dielectric constant is e = 1 + (Nira)/8, where N = 4 is the number of fermion flavors (see 
next section). The inclusion of low order ladder diagrams [llj does not change significantly this value. 

An independent estimate of a in graphene has been obtained from measurements of the carrier-plasmon interaction 
in samples with a finite carrier concentration [26} . The result, a ~ 2.2 is consistent with the previous analysis. 



IV. RENORMALIZATION OF THE FERMI VELOCITY. 1/N EXPANSION 

The first perturbative renormalization results on the running of the coupling constant in graphene were obtained 
in rcf. [T3 on the basis that, even if initially the graphene coupling constant is not small, it runs towards small values 
in the infrared where perturbation theory can be trusted. The logarithmic renormalization described in Sect. Ullis a 
one loop RG result that assumes the coupling constant to be small. 

A standard perturbative framework for strongly coupled fermionic theories is the 1 /N expansion first proposed in 
relation with Quantum Cromodynamics (QCD) [23]. N is the number of fermionic species in the problem that in the 
case of graphene with the spin and valley degeneracies is N = 4. This procedure was followed in the early publications 
[28l [29} and was later retaken in [30I - I32} . 

The simplest non-perturbative calculation amounts to compute the graph in fig. [1] (a) substituting the photon 
propagator by the RPA result which is a resummation of the planar diagrams dominant in the 1/N approximation: 

V(u,\c) = 2 ~ l t2 • ( 12 ) 

From that we get the following equation for the running Fermi velocity [29L [30} : 

1 8vf 8 / arccosp ^4 

where g = (the large N limit amounts to take the limit N —> 00 keeping g fixed). As we see the dependence on 
g is non-perturbative and the growth of the velocity at low energies is slightly different than that in ([5} . 

Another interesting prediction done within this framework is the linear dependence of the quasiparticles lifetime 
with the energy [28} , a distinctive of the marginal Fermi liquid behavior [33} . Experimental evidences of this behavior 
have been described in 0, [H 



V. RENORMALIZATION OF THE FERMI VELOCITY: EXPERIMENTAL CONFIRMATION 



The deduction of the values of the physical parameters from a given experimental measure is one of the most 
important aspects of the renormalization program. This is one of the issues where physics from different areas come 
into play, the best example being the determination of the fine structure constant of QED. By combining solid 
state measurements of the electron precession in magnetic fields or the quantum Hall effect, with theoretical QFT 
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FIG. 2: Experimental measurements of the effective mass of carriers in high mobility graphene samples. Red dots give results 
in suspended samples, and blue dots are for samples on BN substrates. The red and blue lines are fits using eq. 1131 Red: 
eg = I- Blue: eg = 4. The dotted orange line shows a fit obtained neglecting the renormalization of the Fermi velocity. 



calculations of the anomalous magnetic moment of the electron we reach the actual precision of better than one part 
in a trillion [36j | . 

It is fair to say that the Fermi velocity plays in graphene a role similar to the effective mass in the standard Fermi 
liquid theory: it enters in the definition of almost all the phenomenological quantities and, as such, it permeates in the 
value of the experimental observations. Taking for granted a constant, absolute value for this parameter can misguide 
the interpretation of some experiments or hide interesting physical issues. There has been a very recent experimental 
report on the observation of the renormalization of the Fermi velocity [37| confirming earlier more indirect evidences 
This is a very important result which can also affect the interpretation of other experimental observations. 
On the conceptual point of view, the verification of the growth of the Fermi velocity at low energies is interesting 
not only for the condensed matter community but also for the high energy as a beautiful example of renormalization 
at work. Under a more applied point of view, the Fermi velocity being almost the only parameter of the graphene 
model, it affects most of the experimental findings. The recognition that it changes with the energy can offer new 
insights in the interpretation of infrared data. 

The renormalization predicted in eq. [13] is confirmed by the experimental results in 37]. This reference reports 
measurements of the effective mass at different carrier densities for high mobility suspended graphene and graphene 
on BN. These experiments use the temperature dependence of Shubnikov-de Haas oscillations to infer the dependence 
of the Fermi energy, ejr, on the area of the Fermi surface, Sf- The effective mass is defined as 



so that for graphene 



hk F , 

leff = — (15) 
Vf 



A comparison of experimental results an fits based in eq. [13] are shown in Fig. [2] 



VI. OPENING A GAP IN GRAPHENE. EXCITONIC TRANSITION 



Graphene is a semimetal and has zero gap. Under the QFT point of view, the electron mass (gap) is protected by 
the 3D version of chiral symmetry and hence a gap will not open by radiative corrections at any order in perturbation 
theory Quantum electrodynamics (QED) in (2+1) dimensions has been studied at leng th in quantum field theory 
mostly in connection with dynamical mass generation @, I381440T ] and confinement [U, |42| ] Both problems are crucial 
in the physics of graphene, the first is related with the possibility of opening a gap in the system, the second one 
with the screening of the Coulomb interaction. A very interesting review on the QED properties of graphene and its 
possible dynamical mass generation can be found in [43j. 

As mentioned in Sect. UH the coupling constant of QED(3) has dimensions of energy and, in this sense, mass 
generation is greatly facilitated since there is already a mass parameter in the theory. Even then, spontaneous chiral 
symmetry breaking in QED(3) in the 1/N approximation usually requires an unphysical number of fermionic species 



N < 2 [44 1 . 



The situation with graphene is even worse since there is no mass parameter to begin with. The simplest non- 
perturbative calculation that can be done to study the issue of the gap opening in graphene is the RPA type described 
in Sect. IIVI [29|. The absence of a constant term in the inverse electron propagator ensures that no mass is generated in 
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this approximation. The 1/N expansion approach to the problem has been revised recently in [45[ with the conclusion 
that a gap will not open for the physical values of the electronic degeneracy in graphene (N=4). A variational approach 
to the excitonic phase transition in graphene including the renormalization of the Fermi velocity [IgJ also produces 
quite ne gat ive results. A topological analysis of the stability of the Fermi points in graphene against gap opening was 
done in |47| . 

Despite these negative results, many reports on opening gaps in neutral graphene are found in the literature both on 
the experimental and theoretical sides. The possibility of an excitonic insulator was first suggested by Khveshchenko 
48, |49| which also raised the issue of magnetic catalysis [5(| based on earlier studies of QED(2+1) [5l| . 

The possibility of opening a gap when the substrate is commensurate with the lattice was suggested in [U and 
given as a possible explanation of a gap observed in photoemission experiments on epitaxially grown graphene [52j 
and in scanning tunneling microscopy experiments of graphene on graphite [To| . The influence of the substrate in 
this context has also been analyzed in [53|. Monte Carlo calculation in a lattice gauge theory framework predict the 
opening of a gap [U HH ■ A gap has been also claimed in very recent Angular resolved photoemission spectroscopy 
measurements of epitaxial graphene upon dosing with small amounts of atomic hydrogen (55|. Disorder can give rise 
to various strong coupling phases [30L l56l whose physical properties only been partially explored [EH . 

The issue of the gap opening in single layer graphene remains open and some experimental effort is needed to 
dillucidate it. The situation in the bilayer [59|, |6(| is more clear because the gap is not protected there, short range 
interactions can be active at low energies and gaps can open in various ways [60l469| . 



VII. RETARDED COULOMB INTERACTION. EMERGENT LORENTZ COVARIANCE 



As we have discussed, in the static model commonly used in graphene with an instantaneous photon propagator, 
the Fermi velocity grows without bound in the infrared 17] . This sets a lower bound on the validity of the model 



that breaks down when the Fermi velocity approaches the speed of light c. The scale defines an infrared cutoff 5 for 
the theory, and can be computed as: 



fc^exp 



167r(c — v(kn)) 



(16) 



where ItR, vr are the renormalization condition chosen to define the theory A typical estimate will of course produce 
an energy of approximately 100 orders of magnitude below k^. This discussion is very similar to that leading to 
the Landau pole of QED [7(| and simply points to the incompleteness of the theory. In our case there are physical 
low energy bounds much less stringent than that, preventing to access the infrared region. But even if this does not 
impose any real bound on the validity of the model from a experimental point of view, it is interesting to know that 
a complete theory exists of which the static limit is only an approximation. 

In the retarded model pjUZlI the Fermi velocity does not run without stop. The beta function of the Fermi velocity 
has a non-trivial fixed point and the value of vf at the fixed point is precisely the speed of light c. This is a very 
neat example of emergent Lorentz symmetry. 



VIII. DISCUSSION 



The most important parameter in the computation of the observable quantities in graphene is the renormalized 
Fermi velocity that defines, together with the polarizability, the Coulomb coupling constant. 

The program to fix this constant must be similar to the one followed in the case of the electromagnetic fine 
structure constant [36] . In the case of graphene the theoretical determination is much simpler since there are no mass 
parameters relations involved in the calculations. Moreover we have the additional fact that the quantity determining 
Qg (Fermi velocity) is by itself an observable related to the one particle properties of the system. We expect that a 
proper combination of photoemission [H, [34| optical [ll[ and transport [TO, H3, Iz3 measures with the corresponding 
calculations in the QFT renormalization scheme should be enough to determine aa as precisely as needed both 
theoretically and experimentally, similarly to what happens in QED(3+1). 
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